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In the frame of an effective MSSM model without gaugino-mass unification at a grand unification 
. scale, we set a lower bound on the neutralino mass based on the new WMAP data on Qcdm 

' (R-parity conservation is assumed). Our lower bound, m x > 6 GeV, leaves much room for relic 

neutralinos significantly lighter than those commonly considered (m x > 50 GeV). We prove that 
these light neutralinos may produce measurable effects in WIMP direct detection experiments of 

Ch , low energy threshold and large exposure. 
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> ! I. INTRODUCTION 

o 

oo 

■ In a previous paper wc called the attention to relic neutralinos of light masses, i. e. with a mass m x < 45 
GeV. Actually, most analyses on cosmological neutralinos employ a lower bound m x > 50 GeV, i. e. a bound which 
rests on the assumption that gaugino masses are unified at the grand unification (GUT) scale -^QiZT ~ 10 16 GeV. 
However, this unification assumption might not be justified, as already discussed some time ago 013 • Furthermore, 
recent analyses of string models (see, for instance, Ref. |4j) indicate that the initial scale for the running of the 
susy parameters by the renormalization group equations may be much lower than the standard GUT-scale, with the 



I | implication that at the electroweak (EW) scale the gaugino masses may be quite different from what expected in a 
standard SUGRA scheme with GUT-unification assumptions. Thus, for instance, the U(l) and SU(2) gaugino masses, 
c"| Mi and M 2 , may not be related by the standard formula M\ ~ \M 2 at the EW scale. 

It is then quite natural to discuss supcrsymmetric models where M\ and M 2 are considered as independent param- 
eters. Previous papers where schemes of this type are considered include the ones in Refs. HH, HE HSU 113,111 



Eil Eil IbH Eil [la IetI Hit . In Ref. Q we evaluated the neutralino relic abundance VL x h 2 and the neutralino- nucleon 
scalar cross-section cr^^ in an effective Minimal Supersymmetric extension of the Standard Model (effMSSM) 
where the gaugino-mass unification at GUT scale is not assumed. The analysis was performed in a scenario where 
the ratio R = M\/M 2 is smaller than the GUT-unification value Rgut — 0.5, thus allowing the neutralino mass 
m x to be lighter than the commonly- used lower bound of about 50 GeV. We showed that, in the derivation of the 
lower limit on m x , the upper bound on the relic abundance for cold dark matter (CDM), flcDMh 2 , plays a crucial 
role. A lower limit of m x > 5 GeV was derived in Ref. by employing the upper bound flcDiah 2 < 0.3. It was 
also shown how this upper bound on flcDAih 2 is instrumental in providing sizeable values for the neutralino- nucleon 
scalar cross-section at small m x . 

New data on the cosmic microwave background (CMB) 0|i3,El|, also used in combination with other cosmological 
observations, are progressively narrowing down the ranges of the relic abundances for matter (fl m h 2 ) and for some 
of its constituents: neutrinos (fl u h 2 ) and baryons (Qbh 2 ). As a consequence also the range of QcDNih 2 is reaching a 
unprecedented level of accuracy. The importance of an improvement in the determination of SIcdm^ 2 for any cold 
relic particle is twofold. The upper bound {flcDMh 2 ) m ax obviously establishes a strict upper limit for any specific 
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cold species. On the other side, the lower bound (QcDMh 2 ) m in fixes the value of the average abundance below which 
the halo density of a specific cold constituent has to be rescaled as compared to the total CDM halo density. For the 
determination of the rescaling factor £ = p x / po (where p x and po are the local neutralino density and the total local 
dark matter density, respectively), we use the standard rescaling recipe £ = min(l, fl v /t 2 / (r^cr>M^ 2 )min)- In Ref- 19( 
ranges for Q m h 2 and Qbh 2 are derived, by employing CMB data of Refs. |2JJ, |2JJ , 2dFGRS measurements 22( 
and Lyman a forest data j^. From the values quoted in Ref. 01 f° r ^m,h 2 and fl^h 2 and allowing for a 2 — a 
range in UcDAfh 2 , one obtains: {£lcDMh 2 ) m in = 0.095 and (QcDMh 2 ) ma x = 0.131. These are the values we will 
use here. However, one should cautiously still be open to some possible changes in these values as new cosmological 
observational data will accumulate in the future. 

In the present paper we analyse the supersymmetric scenario of Ref. Q in light of these new determinations of 
QcDivih 2 . We first derive approximate analytic expressions which directly relate the lower bound on m x to the upper 
limit on flcDAih 2 (R-parity conservation is assumed). This requires the evaluation of Q, x h? , which is in itself a very 
complicated function of the supersymmetric parameters. In order to extract the leading analytic terms at small values 
of m x , we are guided by a detailed numerical analysis of the supersymmetric parameter space. This will allow us 
to present analytic expressions for the lower bounds on m x which display the link of these bounds to the relevant 
particle-physics and cosmological constraints, such as the lower bounds on sfermions and Higgs-bosons masses and 
the upper limit on flcDAih 2 . Furthermore, on the basis of these results, we show that light neutralinos, with masses 
m x < 45 GeV, may actually be probed by WIMP direct detection experiments with high sensitivities and low energy 
thresholds. 



II. LOWER BOUND ON m x 

The supersymmetric scheme adopted here is the same as the one described in Ref. 0: an effective MSSM scheme 
(effMSSM) at the electroweak scale, with the following independent parameters: M 2 , p, tan/3, rriA, nig, m ; ~, A and 
R = MxjM-2.. Notations are as follows: p, is the Higgs mixing mass parameter, tan/3 the ratio of the two Higgs 
v.e.v.'s, rriA the mass of the CP-odd neutral Higgs boson, mg is a quark soft-mass common to all squarks, mj is 
a slepton soft-mass common to all sleptons, A is a common dimensionless trilinear parameter for the third family, 
At = Ai = Arrig and Af = Amj (the trilinear parameters for the other families being set equal to zero). 

The neutralino is defined as the lowest-mass linear superposition of bino B, wino and of the two higgsino 

states H?, H%: 

X^aiB + a 2 W {3} +a 3 H1 + ai H%. (1) 

Since we are here interested in light neutralinos, we consider values of R lower than its GUT value: Rqut — 0-5; 
for defmiteness, we take R in the range: 0.01 - 0.5. 

We first outline the procedure for deriving analytical bounds on m x from the cosmological upper limit on flcDMh 2 . 
As mentioned in the previous section, the identification of the leading analytic contributions is guided by numerical 
analysis. This is based on a scanning of the supersymmetric parameter space, with the following ranges of the MSSM 
parameters: 1 < tan/3 < 50, 100 GeV < \fj,\, M 2 , m 5 , mr < 1000 GeV, sign(^i) = -1,1, 90 GeV < m A < 1000 GeV, 
—3 < A < 3. The following experimental constraints are imposed: accelerators data on supersymmetric and Higgs 
boson searches (CERN e + e~ collider LEP2 [24| and Collider Detector CDF at Fermilab [l^); measurements of the 
b — » s + 7 decay pif: measurements of the muon anomalous magnetic moment = (g M — 2)/2 |27j (the range 
— 160 < Aa^ • 10 11 < 680 is used here for the deviation of the current experimental world average from the theoretical 
evaluation within the Standard Model; for the derivation see Ref. Q). 

The neutralino relic abundance is given by 



x f 3.3-10- 



n** a = . V/a ~ > (2) 



where < <J ann v > = Xf(a ann v){ nt , (c ann w)i n t being the integral from the present temperature up to the freeze-out 
temperature Tf of the thermally averaged product of the annihilation cross-section times the relative velocity of a pair 
of neutralinos, Xf is defined as Xf = and g+(xf) denotes the relativistic degrees of freedom of the thermodynamic 

bath Ski Xf. For (er a nn v) we will use the standard expansion in S and P waves: (cr ann v) ~ a + -^-b. 
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A lower bound on m x is now derived from Eq. by requiring that 

max ■ (3) 

We first work out approximate analytic expressions for < a ann v > in the regime of small m x {m x < 45 GeV). By 
diagonalizing the usual ncutralino mass matrix in the approximation Mi << A/2, /i, it turns out that light ncutralinos 
have a dominant bino component; a deviation from a pure bino composition is mainly due to a mixture with H^, i. 
e. \ai\ » |o-3 1 >> |ai|, |d4 1 . For the ratio |ci3 S/ l a i i one finds 

f4 ~ sm6 w sin/3 — < 0.42 sin/3, (4) 

|di| w 

where in the last step wc have taken into account the experimental lower bound fx > 100 GeV. 

The dominant terms in (o" ann i>)i n t are the contributions due to Higgs-cxchangc in the s channel and sfermion- 
exchange in the t, u channels of the annihilation process \ + X ~ * / + / (interference terms are neglected). We retain 
only the leading terms in each contribution. Thus, for the Higgs-exchange contribution, dominated by the S-wave 
annihilation into down-type fcrmions, we have, for any final state / — /, 

Higgs „. 2-Ka 2 cr n n n rn 2 f ml [1 — m 2 f /m 2 ] 1 / 2 

(a ann v) f ^a? 9 "^ -Z L Q Q ton 2 /3(l + e/ ) 2 -f * L fl f , (5) 
5 3 sar9 w cos 2 w n% [{2m x y - m\Y 

where Cf is a standard color factor (cj = 3 for quarks. Cf = 1 for leptons), fhf is the fermion running mass evaluated 
at the energy scale 2m x and m/ is the fermion pole mass, e/ is a quantity which enters in the relationship between 
the down-type fermion running mass and the corresponding Yukawa coupling (see Ref. |2sj and references quoted 
therein); in the following evaluations, £/ is negligible, except for the bottom quark, where ef, ~ 0.2. One easily verifies 

- — -~ Higgs 

that when m x < mj, (er ann v) ^ entails a relic abundance exceeding the cosmological bound. 

Notice that < u ann v > turns out to be an increasing function of m x . Then, to obtain a conservative lower bound 
on m x from the condition of Eq. we have first to evaluate an (0, x h 2 ) min which is obtained from Eq. by 
replacing < a ann v > with its maximal value (< cr a „„u >) max > & t fixed m x . In the case of Higgs-exchange contributions, 
this (< o ann v >) m ax is obtained by inserting into Eq. JSJ the maximal value of the product a\a\ tan 2 j3. Taking into 
account Eq. Q and that, for jtia — 90 GeV, the upper bound of tan/3 is 45 [ijj. we obtain (d 2 d| tan 2 (3) max ~ 2.6x 10 2 , 
and in turn 



5-10~ n x f m 2 w 1 [(2m x ) 2 - m 2 A ] 2 

GeV 2 g *( Xf f 2 W (1 + ™ 2 X [1 - m 2 /m 2 ]V2 ' 



(6) 



In the last step of this equation wc have retained only the dominant contribution due to the 6 — 6 final state. 

1/2 

As far as the value of g*(xf) is concerned, we notice that for these light neutralinos Xf ~ 21 — 22, so that 
neutralinos with masses m x ~ 6-7 GeV have a freeze-out temperature Tf ~ Tqcd, where Tqcd is the hadron-quark 
transition temperature of order 300 MeV. For definiteness, we describe here the hadron-quark transition by a step 

1/2 

function: if Tqcd is set equal to 300 MeV, then for m x < 6 GeV one has g*(xf) ~ 4, while for heavier ncutralinos 
g^x f f' 2 ~ 8 - 9. 

The quantity (^ x ^ 2 ) m ^ ss as given by Eq. © is plotted in Fig.l as a function of m x , for the value tua = 90 GeV 
(current experimental lower bound). The solid curve refers to the case Tqcd = 300 MeV; dashed and dot-dashed 
curves denote two different representative values of Tqcd- Tqcd = 100 MeV and Tqcd = 500 MeV, respectively. 
The two horizontal lines denote two representative values for the upper bound on flcDAih 2 : (QcDMh )max = 0.3 
(short-dashed line) and (flcDMh 2 ) m ax = 0.131 (long-dashed line). Fig.l displays how a lower bound on m x is derived 
from an upper limit on flcDAih 2 . In particular, using Tqcd — 300 MeV, one obtains from iyicDMh 2 ) m ax = 0.3 
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the bound m x > 5.2 GcV, a value which increases to m x > 6.2 GeV, when the new value (QcDMh 2 ) ma x = 0.131 is 
employed. Also evident is the variation of the lower bound on m x , when the value of Tqcd is changed. Finally, to 
support the validity of the analytical approximations employed to derive Eq. ijSJ, in Fig.l we also display the scatter 
plot of fl x h 2 , when a numerical scanning of the supersymmetric parameter space is performed. 

We recall that the bound on m x : m x > 6.2 GcV was derived using tua = 90 GeV (which is the present experimental 
lower bound on ttia)- From our previous formulae one obtains that this bound on m x simply scales with itia as follows 



m x [1 - m^r/4 > 5 3 GcV (^^) 2 (7) 

For the sfermion-exchange contribution. (fj an n v) ^ ~ ^y ermlon _|_ _J_^/ e mto j both S-wavc and P-wave 

contributions have to be taken into account. In the regime we consider here, the leading terms due to the down-type 
fcrmions are 



~ sfermion 



a 



8 COS 9 W 



4 m2 A l 

a lCf 



m)/ml 



1/2 



m /,i (2 



f , R ) + 2Y fJ 



lYi,r 
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TO/ 



(8) 
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(9) 



Notations are as follows: Yf t L and are the weak hypercharges for left and right couplings, respectively, mj 
denotes either mj or rriq depending on the nature of the fermion, to ? 1 is the smallest mass eigenvalue for the sfermion 
/, once the mass matrix in the weak-interaction basis /j,, fn is diagonalized. A maximal mixing between the /l, Jr 
fields has been used in deriving Eqs. ||SJ| 0, since these equations are meant to provide a conservative lower bound 
on m x . Notice that for I = e,fi, the Jl ~ Jr mixing is negligible; this case is simply recovered from Eqs. JSJ)-© by 
setting to 



fA 



v) 



sfermion 



With the aid of numerical evaluations, it is found that the leading contributions to (. 

- — - sfermion - — — sfermion - — - sfermion 

J2f ("'aim v)f are provided by the term due to the r lepton: (cr ann v) ~ (<7 ann v) T . This is 

in turn maximized by 



- — - sfermion 
((^anii Vj *)max 



to 2 ./to 2 ,] 1//2 



8 COS 4 9 W 



5 TO T 
2 TOv 



23 
2^7 



(10) 



We denote by (n x h^~ 



sfermion 



the value of Q, x h 2 



derived from Eq. when (cr ann v) is replaced by 

x of Eq. (UJJ. The quantity (^ x h 2 ) s ^" n/lon is plotted in Fig.2, for the value m f = 87 GeV 



(current experimental lower bound); also the scatter plot for the quantity £l x h is displayed, for tua >300 GeV. 



From Fig.2 one finds for m x a lower bound m x > 14 GeV for (f2 



CDM 



0.3, and m x > 22 GeV for 



0.131. The scaling of this last bound with the stau mass is approximately given by 



TO x [l-TO 2 /TO 2 ]V 4 >22GeV (^y 2 - (ID 

In general, one has conservatively to retain as a lower bound to m x the smaller of the two lower limits given 
separately in Eq. J7J) and in Eq. 1|11|) . From these equations one finds that the lower bound of Eq. JJJ is less stringent 
than the one of Eq. (|llf> as long as uia < 2m f . Due to the present experimental bounds on niA and rrif the lower 
absolute bound is the one derived from Eq. J7J), i. e. m x > 6.2 GeV. We parenthetically note that the lower limits 
to x > (15 — 18) GcV found in Refs. are due to the assumption that niA is very large (ttla ~ 1 TcV). 

A plot which shows the transition between the lower bounds established by considering Higgs-cxchange contributions 
and sfermion-exchange contributions is reported in Fig. 3. As was derived by use of the approximate expressions in 
Eqs. 0-tnji the transition point is at m A ~ 200-300 GcV. 
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III. DETECTABILITY OF LIGHT NEUTRALINOS BY WIMP DIRECT MEASUREMENTS 

Now we wish to show that a part of the neutralino population of small m x is indeed cxplorablc with experiments of 
WIMP direct detection with current sensitivities. Let us start by giving an approximate relation between the scalar 
neutralino-nucleon cross-section and the neutralino relic abundance, valid when ^""lar is dominated by /i-exchange 
in the t-channel and a ann by ^-exchange in the annihilation s-channel. From the approximate formulae given in Ref. 
for a"s"aiar° n ^ an d ^x^ 2 wc obtain 

ffi h 2 ) CT ( nuclcon ) ~ i 4 x 10- 40 cm 2 T I ms 2 GeV2 ( UlA 4 fi2) 

l"x'U<W -1-4x10 cm 1 ^ 200 McV J m |[l-m 2 /m|]V2 \m h ) ' ^ 

where is the mass of the lightest CP-even neutral Higgs boson and T is given by 

(03 sina + a4 cosa) 2 (sin a + e s cos (a — j3) sin/3) 2 
(04 cos ,3 - a 3 sin/3) 2 sin 2 /3 (1 + e b ) 2 

In Eqs. H12fl - l|13|) (N\ss\N) is the s-quark density matrix element over the nucleonic state and a is the angle which 
rotates the Higgs fields h[ ^ and into the mass eigenstates h and H . The approximations employed in deriving 
Eas. l|12[l - lll3[l imply a ~ ?, ~ ~ 100 GeV, so that T is of order one. 

Thus, for neutralino configurations with m x < 20 GeV, cr^^ on ^ turns out to be bounded by 

(nucicon) lO^cm 2 GeV 2 

Scalar ~ (ClcDMh 2 ) max m 2 [1 - mf/m 2 ,} 1 / 2 ' 1 ' 

In deriving Eq. O, we have set m s (N\ss\N) = 200 MeV 0. 

The results of a complete numerical evaluations of the quantity ?Cg™ 1 ! i 1 1 con ' ) , where all relevant diagrams for the 
neutralino-nucleon scalar cross-section and for the relic abundance are taken into account, are displayed in Fig. 4. The 
peculiar funnel in the scatter plot for m x < 20 GeV is due to the bound of Eq. (|14f> . 

As was pointed out in Ref. Q, the present upper limits to £ ^scalar ' 1 '' provided by WIMP direct detection 
experiments 

USUI HQ 

do not significantly constrain the supcrsymmctric configurations for the light ncutralinos 
displayed in Fig. 4. Instead, these configurations may be relevant for experiments of direct detection with a low energy 
threshold and a large exposure. An experiment with these features is the DAMA/Nal experiment, whose results, 
after a 4-years running with a total exposure of ~ 58 000 kg • day, show an annual-modulation effect at a 4cr C.L. 
which does not appear to be related to any possible source of systcmatics 0. The analysis carried out by the 
DAMA Collaboration to explain their modulation effect in terms of a WIMP with coherent elastic scattering was 
targeted to a neutralino in the frame of a usual supersymmetric scheme with gaugino-mass unification at GUT, with 
a consequent lower bound on the neutralino mass above 30 GeV. This interpretation was proved to be consistent with 
supersymmetric models with gaugino-unification at GUT |35j| . 

Here we have considered a different supersymmetric scheme which includes significantly lower neutralino masses; 
thus, in order to establish the possible relevance of our low-mass neutralinos for an annual-modulation effect, we 
have to proceed to an extension of previous analyses. To put our arguments into a quantitative basis, we evaluate 
(£ CT scaiar° 0n ^) m ™' dcfmccl as the minimal value of S,^^^ ^ which may produce an annual- modulation effect at n 
standard deviations in a detector with a given exposure S, and for a given velocity distribution function f(v) for relic 
neutralinos in our galaxy, i. e. 



(nucleon)\ n I /-,r\ 

^"scalar > min ~ g (A/) 2 ' 



In Eq. H15|) / is defined as the ratio of the expected direct detection rate, integrated over an energy range (i?!,^), 
to the neutralino-nucleon scalar cross-section ^scalar™''- 
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^ (nucleon) J 
S scalar 

For a monoatomic material of nuclear mass number A, one has 



1 r E ^R(E) dK 



I = N T p 




where Nt is the number of the target nuclei per unit of mass, mjy is the nuclear mass, F(E) the nuclear form 
factor and v m i n is the minimal value of the neutralino velocity to produce an event above the detection threshold. 
Generalization of Eq. (|17l) to non-monoatomic materials is straightforward. In Eq. (|15|l AI is defined as AI = 
(/(June) - /(December))/2. 

We emphasize that, by definition, the condition £^lcaiar — (^scalar* on ly establishes a minimal requirement 

for a neutralino of a given mass m x to be detected in a given experiment. It is obvious that the real detectability of 
the signal depends on the actual amount of the experimental background. For example, an equal amount of signal and 
background would double the value of )min as given in Eq. (|15fl . Estimates of backgrounds are omitted 

here. 

To establish whether our population of light neutralinos may be relevant for the effect measured by the DAMA 
experiment, we have evaluated (Cheater )""« as a function of m x ; the energy range of integration employed here is: 
Ei = 2 keV, E<i = 3 keV, and n is set equal to 4. The nuclear form factors for Na and / are modelled in the Helm form 
[33, with the values of parameters given in Ref. |3t| . In Fig. 4 we give the curves of (£ o , ^"^ 1 ° on ^ ) m m versus m x for a 
sample of different galactic distribution functions (DF), taken among those analysed in Ref. |38|. The intermediate 
curve refers to an isothermal DF with vq = 220 km • s _1 and po = 0.3 GeV • cm -3 (vnjs the local rotational velocity). 
The upper curve refers to a spherical Evans' power-law DF (denoted as A3 in Ref. |3J|) with vq = 170 km • s _1 and 
Po = 0.17 GeV • cm -3 ; the lower curve refers to an axially-symmctric Evans' logarithmic DF with maximal flatness 
(denoted as C2 in Ref. 0]) with v Q = 270 km • s" 1 and p Q = 1.7 GeV ■ cm" 3 . 

From the results displayed in Fig. 4 one sees that indeed a part of our population of relic neutralinos of small m x may 
be relevant for the annual-modulation effect discussed above, especially in case of DF's which entail higher velocities 
and /or higher densities as compared to the standard isothermal distribution. This is a new interesting option which 
adds to the other, still valid, possibility which we discussed in the papers of Ref. .35], on relic neutralinos with masses 
above 50 GeV. 



IV. CONCLUSIONS 



We have considered phenomenological properties of relic neutralinos in a range of low masses (m x < 45 GeV) , which 
is allowed in a MSSM model without gaugino-mass unification at a grand unification scale. We have numerically 
evaluated the relevant relic abundance and the neutralino-nuclcon cross-section, and discussed our numerical results 
in terms of analytic formulae, which display the connections among cosmological properties and particle-physics 
parameters. 

Using the latest determinations of flcDAih 2 and assuming R-parity conservation, we have shown that in a MSSM 
model without gaugino-mass unification the lower bound on the neutralino mass is m x > 6 GeV. We have shown 
analytically how this bound is linked to the cosmological upper bound on flcDAih 2 and to the lower limits on masses 
of Higgs bosons and sfermions. 

The implications of light relic neutralinos with masses < 45 GeV for WIMP direct searches have been analysed. It is 
found that these neutralinos are actually relevant for WIMP direct detection experiments of low energy threshold and 
large exposure. The present results extend to small masses our previous analyses about the effects of relic neutralinos 
with masses above 50 GeV 1351. 
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FIG. 1: Neutralino relic abundance Q x h 2 as a function of the mass m x . The solid curve denotes (fix^lmi™ 8 given by Eq. © 
for Tqcd = 300 MeV. Dashed and dot-dashed curves refer to the representative values Tqcd = 100 MeV, Tqcd = 500 MeV, 
respectively. The two horizontal lines denote two representative values of QcDiuh 2 : Q,cdmIi 2 = 0.3 (short-dashed line) and 
QcDNih 2 = 0.131 (long-dashed line). The scatter plot is obtained by a full scanning of the supersymmetric parameter space. 
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FIG. 2: Neutralino relic abundance Q, x h as a function of the mass m x . The solid curve denotes {n x h 2 ) s J^ mion derived from 

^ sf ermion 

Eq. (|2j when (o- a nn «) is given by the maximal value of (cr ann v) (see Eq. CDJ). Tqcd is set equal to 300 MeV. The 

two horizontal lines denote two representative values of QcDAih 2 : QcDMh 2 = 0.3 (short-dashed line) and QcDAih 2 = 0.131 
(long-dashed line). The scatter plot is obtained by a full scanning of the supersymmetric parameter space with ttia > 300 GeV. 



10 



A. Bottino, PI Dtmato, K: Fornenfio, 3, Scopel (SO 03) 




4 5 6 7 8 910 20 30 40 50 



n\ (GeV) 



FIG. 3: Dashed and solid curves give the variation of the lower bound on m x as a function of tua for {p.cDMh 2 ) max = 0.3 
and for (QcDAih 2 ) ma x = 0.131, respectively. For each value of (0.cDi\ih 2 ) ma x the region on the left of the relevant curve is 
forbidden. The two lines are derived from the analytical expressions obtained in the text for the contributions to (cr ann v) due 
to Higgs-exchange and sfermion-exchange. Also displayed is the scatter plot of a full numerical scanning. 
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FIG. 4: Scatter plot of £ CT icaiar ° n ' versus m x . Crosses (red) and dots (blue) denote neutralino configurations with Q x h 2 > 
(yicDMh 2 ) m in and Cl x h 2 < (flcDMh 2 ) m i n , respectively. The curves give the sensitivity line, (£cr£^° °') m j„ versus m x , for 
a Nal detector, whose features are discussed in the text. The intermediate curve refers to an isothermal DF with vo = 220 
km ■ s _1 and po = 0.3 GeV ■ cm -3 . The upper curve refers to a spherical Evans' power-law DF (denoted as A3 in Ref. plijo 
with vo = 170 km ■ s -1 and po = 0.17 GeV • cm -3 ; the lower curve refers to an axially-symmetric Evans' logarithmic DF with 
maximal flatness (denoted as C2 in Ref. [jjl ') with vo = 270 km • s _1 and po = 1.7 GeV ■ cm -3 . 



